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Abstract: Proton-proton fusion is considered in the effective field theory of Kaplan, Savage 
and Wise. Coulomb effects are included systematically in a non-perturbative way. Including the 
dimension-eight derivative coupling which determines the effective ranges of the scattering ampli- 
tudes, next-to-leading order corrections to the fusion rate are calculated. When the renormalization 
mass is well above the characteristic energies of the system, this contribution gives a rate which 
is eight percent below the standard value. The difference can be due to an unknown counterterm 
which comes in at this order. 

During the last year much progress has been made in understanding the low-energy 
properties of few-nucleon systems from the effective field theory of Kaplan, Savage and 
Wise[||]. Both elastic and inelastic processes can be considered. Bound states like the 
deuteron can be treated via its interpolating field which replaces the need for explicit 
wavefunctions^]. A more general review of the whole approach has recently been given 
by Kaplan @. 

Including the two lowest order interactions with coupling constants Co and C2, the 
effective Lagrangian for the nucleon field N T = (p, n) with mass M can be written as 

C = iV+ (idt + N - C (N T IIN) • (iV T niV)t 

+ ^C 2 {{N T V 2 UN)-{N T UN^ + h.c.} (1) 

where the operator V = ( V — v)/2. The projection operators LTj enforce the correct 
spin and isospin quantum numbers in the channels under investigation. More specif- 
ically, for spin-singlet interactions 11, = cr 2 T2Tj/ 1 \f% while for spin-triplet interactions 
Hi = 0"20"iT2/\/8- This theory is now valid below an upper momentum A which will be 
the physical cutoff when the theory is regularized that way. For momenta much smaller 
than the pion mass, we can consider the pion field integrated out and all its effects soaked 
up in the two coupling constants Co and C2. Then the value of the cutoff A will be set 
by the pion mass m n . In this momentum range all the main properties of few-nucleon 
systems should then in principle be given by the above Lagrangian. More accurate results 
will follow from higher order operators in this field-theoretic description 0]. 

1 On leave of absence from Institute of Physics, University of Oslo, N-0316 Oslo, Norway 



1 



From calculations of proton-neutron elastic scattering using the PDS regularization 
scheme]!]] [§] or the equivalent OS scheme ||, the a priori unknown coupling constants Co 
and C2 can be determined in terms of experimental quantities. In the spin-triplet channel 
the deuteron will appear as a bound state and the corresponding lowest order renormalized 
coupling constant is 

c »<"> = f,^ < 2 > 

where 7 = 45.7 MeV is the momentum in the deuteron bound stateQ 01. The renormal- 
ization mass p can be chosen freely in the interval 7 < p < m n and physical results should 
be independent of its precise value. Similarly, for the dimension-eight coupling constant 
one finds 

^rf/.A 47r ( 1 \ 2 Pd 



where pd = 1.76 fm is the spin-triplet pn effective range scattering parameter evaluated at 
the deuteron pole. Together with the similar coupling constants in the spin-singlet channel, 
many properties of the proton-neutron system can be calculated. In particular, results for 
proton-neutron radiative capture n + p — > d + 7 have recently been obtained^] || which 
from the hadronic point of view is very similar to proton-proton fusion p + p — * d + e + + v e . 

The coupling constants in the pp channel can also be matched to experimental scat- 
tering data at low energies. However, in this energy range the Coulomb effects become 
important and must be separated out. This has recently been done within the frame- 
work of the same effective theory fjj and also within standard quantum mechanics ||. The 
leading order coupling constant in the spin-singlet pp channel can now be written as 

= %ttt\ — ^ 

M l/a(p) — [i 

where a would be the strong scattering length when there were no Coulomb interactions. 
But these turn it into a //-dependent quantity which can be determined from the measured 
scattering length a p = —7.82 fm using 



1 1 



a(/j.) a p 



ln ^ + 1 _3 

aM 2 



(5) 



where a is the fine-structure constant and Ce = 0.5772 ... is Euler's constant. Since 
the proton mass is so much heavier than l/a p , the Coulomb correction is seen to be 
surprisingly large. This has been known for a long time and was previously expressed by 
the corresponding Jackson-Blatt relation of the same form and obtained from potential 
models [|l(]]. The higher order coupling constant C\ in this channel will have the same 
form as C2 in (||) but with 7 replaced with l/a(p) and with the proton-proton effective 
range p p = 2.79 fm. We have shown that it is not affected by Coulomb corrections to this 
order in the effective theory. However, the C\ coupling gives an important contribution to 
the scattering length (||) which picks up an additional term —pp p /2 in the parenthesis]^]. 
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Including only the leading order couplings Cq'^, we have recently calculated the rate 
for proton-proton fusion p + p — > d + e + + v e as it takes place in the Sun||. Since the 
initial energy E = p 2 /M is then so low, it can be taken to zero. Thus also the momenta of 
the final state leptons can be ignored. Stating the result in terms of the standard reduced 
matrix element A(E) which is dimensionless [11] ]i~2| , we obtained the result 



A (0) = e x - 2aMa p I{x) (6) 
where the parameter x = aM/7. Here we have introduced the function 

I(x) = --e x E 1 ( X ) (7) 

where E\ (x) is the exponential integral function. This is in full agreement with the corre- 
sponding result obtained in the zero-range approximation of nuclear potential models [fig]. 
With the above values for the different parameters, we have x = 0.15 and 1(0.15) = 4.96. 
The matrix element then becomes Ao(0) = 2.51 so that in the rate we have Ag(0) = 6.30. 
This is 10% below the value A^ ot (0) = 7.0±0.05 obtained in the most complete calculations 
including higher order effects [fl3|[ E 



When we now include the effects of the next order coupling C2 to the fusion process, it 
gives rise to the Feynman diagrams shown in Fig.l. In the chain of bubbles connected by Co 
interactions, the two protons interact via the Coulomb potential Vc(t) = a/r. In a single 
bubble the particles propagate from zero separation and back to zero separation. The 
contribution of a single bubble in the chains is thus Jo(p) = Gc(E; 0, 0) where Gc{E; r, r') 
is the Coulomb propagator in coordinate space of the two protons with energy E. We have 
found it most convenient to calculate the value of the bubble quantity Jo(p) in momentum 
space. With J = fd 3 q/(2ir) 3 it is 

J (p)= / / G c (E;q,q!) (8) 

Jq Jq' 

when expressed in terms of the Fourier-transformed Coulomb propagator Gc(E; q, q'). 
This integral can now be done using either dimensional regularization in the PDS scheme 
or a standard momentum cutoff to regularize the ultraviolet divergence it contains(?j[||]. 
The result has been used to obtain the Coulomb-corrected scattering length (|5]) and the 
leading order fusion result (||). 

In some of the diagrams in Fig.l the derivative operator V acts on such Coulomb- 
dressed bubbles and gives rise to more divergent integrals. One of the simplest is 

J 2 (p)= [ [ q 2 Gc(£;q,q') (9) 



Using a method based upon the the Fourier-transformed Coulomb wavef unctions, we have 
obtained a finite result for this integral|7|. However, the method lacks a firm foundation 
and cannot easily be applied to other integrals of a more general form. We will therefore 
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Figure 1: Feynman diagrams contributing to the fusion process to first order in the derivative 
coupling Ci denoted by black square. The wiggly line indicates the weak current while the cross 
circle denotes the action of the deuteron interpolating field. 

here use a more direct and simpler method based upon the functional equation satisfied 
by the Coulomb Green's functionFJ. Introducing the free propagator 

G (£;q,q') = 2 M „ . (27r) 3 £(q-q') (10) 

the functional equation is Gc = Go + GqVcGc- The integral @ can then be written as 

.72 (p) = P 2 Mp)+ I I (q 2 -p 2 )G (£7;q,q / ) 

+ 1111 (^-V 2 )G {E;^)Vc{^)G c {E^M) 

Jq Jq' Jk A' 

2 This new method arose out of a discussion with Peter Lepage to whom we are much thankful. 
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where Vc*(k, k') = 47ra/(k — k') 2 is the Fourier transform of the Coulomb potential. The 
first integral is now zero by dimensional regularization. Integrating over k in the second, 
we see that the first two factors just give —M. In the integration of Vc(q, k') over q, we 
can now shift the integration variable q — > q — k' and use 

Aira 

7T = OtfJL (11) 

q q 2 

This integral is also zero with dimensional regularization, but it contains a PDS pole in 
d = 2 space dimensions which gives the non-zero contribution. The remaining integrations 
over q' and k' give then finally simply Jo(p). Collecting all the factors, we thus obtain the 
simple result 

Ji(p) = (P 2 ~ a(iM)J (p) (12) 

This agrees with the result using our previous method except for terms which are smaller 
by factors of the order of aM/ \x. These can in practice be neglected since we always take 
the renormalization point [i 3> aM. 

Using the same method to evaluate other similar integrals arising from the Feynman 
diagrams in Fig.l, we can then simplify and gather all the contributions into a modified 
reduced matrix element. A part of this new contribution goes into the next-to- leading 
order scattering length a p in the lowest order result (^). The remaining terms simplify to 



M'y 2 C p -4- C d 

A>10) C 2 rf ^( / x- 7 )Ao(0)-a p7 2 (/i-7)^r A 



However, there is also a wavefunction renormalization constant yff which to this order 
will modify the lowest order result @. It will enter in all calculations involving the bound 
state deuteron and has been calculated by Kaplan, Savage and Wise||. In this version of 
the effective theory without pions it is given by 

X> = 1-C^^-J)(V-2J) (14) 
The full result for the reduced matrix element in next-to-leading order is thus 

A NLO = ^L+A 2 (15) 

Expanding this now to first order in C| with the renormalized value (|3|), we get the final 
result 

Anlo(0) = Aq(0) (l + - 1Pd j - a pl \^ - 7)^^ (16) 

The last term is seen to be dependent in general on the renormalization mass [x. From a 
physical point of view, the result should be independent of this arbitrary parameter. What 
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makes this possible here, is the presence of a new, local interaction which comes in as a 
counterterm at this order of perturbation theory. It will be /^-dependent in such a way as 
to make the overall result independent of p. A very similar situation arises in the process 
n + p — ► d + 7 where such a counterterm also is present § . The a priori magnitudes 
of these counterterms are determined by physics on scales shorter than included in the 
effective theory. An absolute prediction of the proton-proton fusion rate is thus not possible 
at this next-to-leading order as long as this counterterm is not determined by other means. 



Here we will instead compare our result (16) with the corresponding result from po- 
tential models. When we take /i > 7, the dependence of the result on this arbitrary 
renormalization mass becomes negligible and we find 

Ajvlo(0W> 7 = A (0) + ^ 7 Pd) + \a vl 2 ( Pv + p d ) (17) 

With the previous values of the different parameters, we obtain for the reduced matrix 
element the value A(0) = 2.54 which is just a 1.4% addition to the leading order result. 
This is surprisingly small, but results from an almost total cancellation between the two 
effective-range corrections in (17). In the full rate, it corresponds to a value which is 8% 



below the accepted value from nuclear potential models] 13] |14|]. However, the structure 
of our result is very similar when compared to what one obtains in the corresponding 
effective range approximation [TJ]. The last term in (|l7]) is then exactly the same, while in 
the first term the factor (1 + \^pd) is replaced by (1 — 7Pd)~ 1 / 2 which in nuclear models 
is the normalization factor of the deuteron wavefunction. To lowest order in the deuteron 
effective range parameter p d this is then just the same. But 7/^ = 0.41 is really not a small 
expansion parameter and higher order terms in the expansion of the square root give a 
sizeable contribution. These can only be reproduced in the effective theory when including 
even higher order interactions. At the same time this will then also bring in new and 
unknown counterterms. It thus seems difficult for this effective theory to compete with 
the more accurate results obtained from potential models where the deuteron effective 
range corrections are included to all orders. 

A more detailed presentation of the calculation and discussion of our results including 
the presence of the counterterm, will be presented elsewhere. 

We want to thank the organizers of the workshop on "Nuclear Physics with Effective 
Field Theory" at the Institute of Nuclear Theory during which this work came to conclu- 
sion. In addition, we are grateful to the Department of Physics and the INT for generous 
support and hospitality. Xinwei Kong is supported by the Research Council of Norway. 
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